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Abstract 

High redshift radio galaxies (HzRGs) are key targets for studies of the formation 
and evolution of massive galaxies. The role of dust in these processes is uncertain. 
We have therefore observed the dust continuum emission from a sample of z > 3 
radio galaxies with the SCUBA bolometer array. We confirm and strengthen the 
result found by Archibald et al. (1), that HzRGs are massive starforming systems 
and that submillimeter detection rate appears to be primarily a strong function of 
redshift. We also observed HzRG-candidates which have sofar eluded spectroscopic 
redshift determination. Four of these have been detected, and provide evidence that 
they may be extremely obscured radio galaxies, possibly in an early stage of their 
evolution. 
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1 Introduction 



The discovery of a tight relation between the masses of the stellar bulges of 
galaxies and the masses of their central black holes (0) has reemphasized the 
importance of high redshift active galactic nuclei (AGNs) as probes of galaxy 
formation. They are therefore the target of many studies of starformation in 
their host galaxies. Since a large proportion of the optical and UV emission 
from a starburst is absorbed and reprocessed by dust to restframe far-IR 
wavelengths, a good way to investigate the roles of dust and star formation in 
these distant galaxies is to measure their (sub-)millimeter continuum emission. 
Another advantage of observing at those wavelengths is that dimming due to 
cosmological distance is compensated for as the peak in the dust emission is 
redshifted into the observed passband (e.g.,0) 
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The formation of massive elliptical galaxies may be best studied by tracing the 
evolution of powerful radio galaxies. They require the presence of an actively- 
accreting super massive black hole (~ 10 9 M Q ; 0), are very extended and 
luminous at many wavelenghts allowing detailed studies, and their selection 
at radio wavelengths is not affected by the presence of dust. 



Archibald et al. (1) have conducted a submillimeter survey of 47 radio galaxies 
over an interval of 0.7 < z < 4.4. They found a substantial increase in 850 /mi 
detection rate with redshift and that the average 850 /zm luminosity L 850 rises 
at a rate (1 + z) 3-4 out to z ~ 4. While it is true that a buried AGN could be 
responsible for heating the dust, there are strong indications that starbursts 
(SFR ~ 1000 M yr _1 ) are the dominating mechanism: (i) there is no strong 
correlation between radio power and submillimeter luminosity (QUI), (ii) for 
all hyperluminous infrared galaxies the starburst outshines the AGN at rest- 
frame wavelengths longwards of 50 //m (0;0)) an d (in) the extent (up to 30 kpc) 
of the CO and dust-emission found in 3 HzRGs (j8j |9J) is not readily explained 
by a central process. 

Until recently only few HzRGs at z > 3 were known, hampering statisti- 
cal studies. However, using radio and near-IR selection we have significantly 
increased the sample of z > 3 HzRGs at Keck (e.g., fToh ) . Here we present 
submillimeter observations of these. Additionally, we have observed an in- 
triguing class of HzRG-candidates for which no spectroscopic redshift could 
be determined. One of these candidates had been studied before in detail 



(WN J0305+3525: llll) and was found to be very submillimeter luminous. 



2 Sample Selection and Observations 



We have selected all known HzRGs with redshifts z > 3 without prior SCUBA 
observations. We aimed to observe a significant sample of HzRGs to comple- 
ment the observations of Archibald et al. (|l|) and obtain better statistical info 
at the highest redshifts. The sample is the result of an ongoing effort by our 
group and others (e.g..lloT) to find distant radio galaxies based on Ultra Steep 
Radio Spectrum (USS) and near-IR identification selection criteria. Also 10 
HzRG-candidates have been observed that did not show optical emission lines 
in deep spectroscopic observations at Keck. 

The observations have been carried out over 4 semesters in 4 years (August 
1999 - January 2002) with SCUBA £J at the James Clerk Maxwell Tele- 
scope. We observed at 450 /im and 850 /im wavelengths in photometry mode. 
The optical depth r§5o varied between 0.14 and 0.38 with an average value 
of 0.26. The data were reduced following standard procedures outlined in the 
SCUBA Photometry Cookbook. 
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3 Results and Discussion 



3.1 HzRGs with z > 3 



In the course of our program we observed 23 HzRGs, of which 12 were detected 
with a signal-to-noise ratio larger than 3. Adding these data to the sample of 
Archibald et al. (0) increases the number of 3a-detections at z > 3 from 5 to 
16. There was no attempt to correct for a contribution to the submillimeter 
continuum from the non-thermal radio emission, since most sources were se- 
lected using an USS criterion. Moreover, extrapolation from the radio regime 
to submillimeter wavelengths is very uncertain at best, and is likely to result 
in an overestimate of the non-thermal flux density due to steepening of the 
radio spectrum. For this reason (and following UM we have chosen to use the 
"un-corrected" values from Archibald et al. (0). Figure 1 shows two histograms 
representing the number of HzRGs that were observed at each redshift and 
the number of HzRGs that were actually detected given a 2a or 3a detection 
criterion. Regardless of the actual detection limit used, a marked increase in 
the fractional detection rate from ~ 15% at z < 2.5 to ~ 50% at z > 2.5 is 
apparent, confirming the result found by Archibald et al. (|l|). 



We have correlated the inferred submillimeter luminosities with radio power 
to search for signs of significant heating by the AGN. No strong correlation 
is apparent (Fig. 2). Recently, it was found that the size of the radio source 
and L 850 are correlated for a sample of radio-loud QSOs at z = 1.5, with 
smaller (i.e. younger; Eil ) sources being more luminous (14). This effect was 
attributed to a possible relation between the jet-triggering event and a short- 
lived starburst or quasar-heated dust in QSOs. The effect was not found for a 
matched sample of radio galaxies at z ~ 1.5. We have checked whether there 
is a correlation between Lg5o and radio source size (age) in our larger sample. 
Figure 2 shows that they are not strongly correlated. However, size of the 
radio source depends not only on age but also on jet-power and enviroment 
which is hard to correct for in our heterogeneous sample. Therefore we cannot 
exclude that L 850 and age of the radio source are correlated. 



The (sub-) millimeter emission of starbursts should be largely indepent of view- 
ing angle. Unification schemes (e.g., EH) therefore suggest that we can combine 
observations of HzRGs with those of QSOs, allowing us to study the evolu- 
tion of starformation in their host galaxies over a larger redshift range. Figure 
3 shows the average observed (sub-) millimeter fluxdensitv as a function of 
redshift for mostly radio-quiet QSO samples |l6l Il7; lit lil ) and HzRGs. Re- 
assuringly, the trends mimic each other and they indicate that the average 
observed flux density rises to redshifts of z ~ 4 after which it declines. 
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Fig. 1. Histograms of 850 /mi SCUBA detections (dashed) with S/N > 2 (left) and 
S/N > 3 (right) versus the total number of observed radio galaxies (grey) as a 
function of redshift. The data are from with the detections from our program 

shown in black. At z > 2.5 ~ 50% — 67% of the galaxies are detected, as opposed 
to ~ 15% at z < 2.5. 
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Fig. 2. The left plot shows Lsso versus radio power. The right plot shows Lgso versus 
linear size. No obvious correlations are present. 



3.2 Obscured "no-z" HzRGs 



During our search for HzRGs among USS sources, we found that 33% of the 
selected canditates were at redshifts z > 3 (10). However we also found that 
~ 24%, mostly with a compact radio morphology, do not show clear signs of 
emission lines and that a significant fraction (~ 10% of the entire sample) 
does not even show continuum emission. Possible explanations for this may 
be that (i) these are galaxies in the "redshift desert" (1.5 < z < 2.3) where 
only weak emission lines fall in the optical window, (ii) they are pulsars, (in) 
they are extremely obscured AGN (iv) or at such high redshift (z > 7) that 
Lya is shifted out of the optical passband. 
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Fig. 3. The average measured (sub-)millimeter fluxdensity binned in redshift has 
been plotted for QSOs and radio galaxies. The squares represent observations of 
radio galaxies at 850 /mi, the triangles QSOs at 850 /mi, while the circles represent 
observations of QSOs at 1250 fim multiplied by a factor of 2.5 for easier comparison. 

Figure 4 shows one of these objects, WN J0305+3525, which was studied in 
detail (jilt ). It is a strong (sub-) millimeter detection with S 850 = 12.5 ± 1.5 mJy 
and S1250 = 4.2 ± 0.6 mJy, and appears associated with a number of faint 
K ~ 21 mag and J > 23 mag objects. The radio morphology is compact (9 < 
1.9 "), as is typical of the other candidates. Based on circumstantial evidence 
like the near-IR magnitudes and morphology, the (sub-)millimeter ratios and 
compact radio source morphology it was argued that WN J0305+3525 is a 
highly obscured HzRG possibly at a redshift z = 3 ± 1 in a young stage of its 
evolution ffllh . We have detected 3 more of these HzRG-candidates with flux 
densities of S$5o ~ 5 ± 1 mJy resulting in a detection rate of 4 out of 10, very 
similar to what is found for the confirmed HzRGs. This further supports the 
idea that they may be young, extremely obscured HzRGs. 



4 Summary 

• Submillimeter detection rate appears to be primarily a function of redshift. 
If this is interpreted as being due to a change in the intrinsic far-IR lumi- 
nosity, it would be consistent with a scenario in which the bulk of the stellar 
population of radio galaxies forms rapidly around redshifts of z = 3 — 5 after 
which they are more passively evolving. 

• A significant fraction (10%) of USS sources may be very obscured radio 
galaxies. The obscuration may be related to their young evolutionary stage. 

• In contrast to the submillimeter sources, HzRGs have accurately determined 
redshifts and host identifications. Even with the advent of SIRTF this is 
not likely to change. Therefore, HzRGs are a key population for studies 
of galaxy formation in the early universe, allowing detailed follow-up mm- 
interferometry observations to study their dust and gas content. 



5 




3 h 5 m 47.60 s 47.40 s 47.20 s 47.00 s 3 h 5 m 47 60 s 47 40 s 47.20 s 47.00 s 

Right Ascension (J2000) Right Ascension (J2000) 



Fig. 4. Keck/NIRC K-band (left) and J-band (right) images with 4.85 GHz VLA 
radio contours overlaid of the HzRG candidate WN J0305+3525. Note the multiple 
components; objects 2 and 4 are not detected in the 58min J-band image. The 
bright object to the SW is a spectroscopically confirmed star. The circle with 3" 
radius represents the nominal 3<r astrometric uncertainty for the centroid of the 
850 /im emission and dashed lines indicate the LRIS slit positions. 

The work of M.R. and W.v.B. was performed under the auspices of the U.S. 
Department of Energy, National Nuclear Security Administration by the Uni- 
versity of California, Lawrence Livermore National Laboratory under contract 
No. W-7405-Eng-48. 
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